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Abstract. It is shown that the probability of quantum- mechanical transmission 
across a phase space bottleneck can be compactly approximated using an operator 
derived from a complex Poincare return map. This result uniformly incorporates 
tunnelling effects with classically-allowed transmission and generalises a result 
previously derived for a classically small region of phase space. 
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1. Introduction 

There has recently been a resurgence of interest in tlie classical transition state theory of 
molecular reactions. Results that were historically restricted to two degrees of freedom 
m E] have been generalised to arbitrary dimensions using the construction of normally 
hyperbolic invariant manifolds (or NHIM's) [3]-[7j. It is natural to ask how classical 
structure such as the NHIM is reflected in the quantum-mechanical problem, which 
corresponds to scattering from a multidimensional potential barrier. 

An answer to this question has been offered in jH], where a description is given 
of quantum mechanical transport across a phase-space bottleneck using dynamics 
linearised around a certain complex periodic orbit. In using linearised dynamics these 
results are restricted to a classically small region of phase space and energies that are 
no larger than 0{h) above a transmission threshold. In this paper it is shown how fully 
nonlinear dynamics may be incorporated in this approach, resulting in a description 
of transport which is not restricted to a small region of phase space or energy range 
above threshold. The current approach is based on a quantisation of a classical normal 
form Hamiltonian, although the final form can be expressed in such a way that explicit 
calculation of a normal form is not necessary. 

To describe the result more concretely, let us consider a waveguide problem 
with configuration space coordinates (x, y) in which x represents longitudinal position 
along the waveguide and y represents transverse vibrations. If necessary we can let 
y = {yi, ■ ■ ■ ,yd) be multidimensional. In chemical applications, x might be a reaction 
coordinate with x large and negative corresponding to decoupled reactant molecules and 
X large and positive corresponding to decoupled product molecules (see Figure Hj), while 
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Figure 1. A schematic illustration is given of the waveguide problem considered in 
this paper, showing the case for an energy above threshold in which the transition 
state (part of which is indicated schematically by a dashed curve) separates reactants 
from products where the waveguide is at its narrowest. The surface shown is meant 
to illustrate a level surface in configuration space of the potential energy in the case 
of a three-degree-of-freedom Hamiltonian that is of the form kinetic-plus-potential 
(although no such explicit assumption regarding the form of the Hamiltonian is made 
in the calculation). Note that we use the term "waveguide" in a generalised sense in 
this paper to mean a Hamiltonian with a phase space bottleneck which asymptotes to 
a vibrational problem decoupled from free motion, and we do not deal explicitly with 
hard-wall or discontinuous potentials. 



y describes internal vibrations of the reacting molecules. The quantum mechanics of 
this problem are described using a scattering matrix, which we write in the form 




where, for example, the block tp^ maps asymptotic incoming states on the reacting 
side to the corresponding asymptotic outgoing states on the product side. In chemical 
jargon, tpR gives state-selected reaction rates (labelled by the incoming mode number) 
together with the distribution of product states (labelled by the outgoing mode number). 
In this paper we will describe a semiclassical approximation for the operator 

n{E) = t^p^tpR 

which gives a probability of transmission for states incoming on the reactant side but 
which sums over outgoing states and does not give the distribution of product states. 
Although containing less information than the scattering matrix, this reaction operator 
has clear experimental relevance and, importantly in this context, admits semiclassical 
approximations which are considerably simpler. As described in detail in |H], this is 
because the orbits used in semiclassical approximation of S{E) are singular near the 
boundary of the reacting subset of phase space whereas those used for IZi^E) are not. 

The reaction operator 'R-{E) has a clear relationship with the geometry of the 
classical transition state. It acts on the Hilbert space of asymptotically propagating 
incoming states and the classical analogue of this space is a Poincare section obtained 
by fixing the reaction coordinate and the total energy E, for which {y,Py) provide 
canonical coordinates. Let E be greater than threshold so that there is a nonempty 
reacting subset V of — the boundary of V is the intersection with S"^ of the stable 
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manifold of the NHIM. Then a phase space representation of 7^(£') such as the Weyl 
symbol W-j^{y,Py) tends to the characteristic function of the reacting region V 

^niy,Py) ~ Xv{y,Py) 

in the classical limit. Moreover, for finite values of h, 'W^{y,py) also incorporates 
quantum effects such as tunnelling, especially important outside V and near its 
boundary. 

An explicit semiclassical approximation was presented for iZ{E) in [Hj, of the form 

nE) = (1) 

where the operator ^{E) is constructed from the linear stability properties of a complex 
periodic orbit. This complex periodic orbit has a real initial condition in the interior 
of V and returns to it after encircling the transition state region in a net imaginary 
time. The formula was derived by using a separable approximation of the Hamiltonian 
in the transition-state region to match waves propagating in the reactant and product 
channels. This separable approximation is valid only insofar as the transition state is 
small on classical scales and the result should therefore work only when the energy E 
is within 0{h) of threshold, where the Liouville volume of V is 0{h'^). 

We will now show that, as long as the operator ^{E) is interpreted using fully 
nonlinear dynamics in a neighbourhood of the periodic orbit. Equation (pj is in fact 
valid over classical scales. The difference between the approach in [Hj and the philosophy 
applied here has an analogy in the classical treatment of one-dimensional WKB solutions 
near turning points. The simplest way to treat turning points is to approximate the 
potential using a truncated Taylor expansion (linear for a single turning point and 
quadratic for two coalescing turning points) and to use the resulting solutions to match 
standard WKB approximations on either side. The method of comparison equations 
in], on the other hand, seeks a change of variable which (up to higher-order corrections 
in h) maps the potential more globally into a linear or quadratic form as required and 
this has the advantage of giving uniform results which apply over classical length scales. 
The approach in jH] is analogous to the method of truncating Taylor series whereas in 
this publication we pursue a transformation into normal form that is similar in spirit to 
the method of comparison equations. 

The difference is that, for multidimensional problems, a deformation of 
configuration space variables alone as used in the method of comparison equations 
does not have sufficient range to put the problem in a solvable form and we must 
use transformations in phase space fOl El 1^- In this way a Hamiltonian which 
is a quantum analogue of the normal forms in [4| can be used and we arrive at a 
problem which, while not separable, is simple enough that scattering solutions can be 
written down. More importantly, the information we need to approximate 'R-{E) can be 
formulated in such a way that explicit reference to the normal form can be removed and 
the end result is a formula (of the same form as Equation (P)) which can be understood 
simply in terms of complex orbits starting and finishing on T}^. 
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We conclude this section with a brief overview of the paper. The essential features 
of the classical normal form are described in Section |21 and an overview is given of the 
corresponding quantum Hamiltonian. The normal form Hamiltonian is not separable 
but does have scattering solutions that are of separable form and these are described 
in Section El Because these scattering solutions do not fully separate in the eigenvalue 
equation and because the normal form Hamiltonian is not of kinetic-plus-potential type, 
using them to describe the transmission properties of a general scattering state is not 
straightforward. Nevertheless a simple solution to this problem is possible, which is 
outlined in Section HJ The final step in obtaining a usable result is to present the 
scattering solution in a basis-independent way, which we do in Section El in terms of 
quantised complex Poincare mappings. Conclusions are presented in Section IHl 

2. Classical and quantum normal forms 

The basis for the calculation in this paper is a quantisation of the classical normal 
form for the Hamiltonian around an equilibrium. The classical normal form and its 
connection with classical transition state theory are described in detail in In this 
section we will describe the essential results and adapt some of the notation for our 
own purposes. We will then describe the important properties of a quantisation of this 
normal form. 

2.1. The classical normal form 

Let canonical coordinates {qo,Po, q,p) = {qo,Po, Qi, ■ ■ ■ , Id, Pi, ■ ■ ■ ,Pd), be chosen so that 
the quadratic part of the Hamiltonian is 



We denote by / = 1 + d the number of degrees of freedom and we will refer to 
(go,Po) as the reaction coordinates. In the context of coUinear molecular collisions 
it is useful to let go = define a dividing surface between reactants and products with 
go < corresponding to reactants and go > corresponding to products. The reaction 
coordinates (go,Po) are useful in interpreting the dynamics in this system, in which 
reaction amounts to crossing a parabolic potential barrier in that degree of freedom. It 
will also be useful for computational purposes however to allow alternative coordinates 
(g, P) defined rotation of (go,Po) so that 



We will refer to / as the reaction action. It is positive for nonreacting trajectories, which 
are repelled by the parabolic barrier, and it is negative for the reacting trajectories, 
which cross over (see Figure |21). 




i=l 



I=\{ql-pl)=QP. 
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Figure 2. In classical phase space, we can identify in and out rcactant and product 
channels with sectors in the Q-P-plane as illustrated in this figure. Orientations here are 
for dH/dl < 0. The dashed line labelled R/P divides reactants from products and the 
dashed line labelled in/out divides incoming from outgoing trajectories. For example, 
a point in phase space which projects to the sector in the QP-planc between dashed 
lines and containing the positive Q-axis is in the incoming reactant channel. The 
hyperbolae illustrate (projections of) typical trajectories evolving from the incoming 
to the outgoing channels. Note that in the classical picture dynamics is limited to an 
energetically allowed region bounded by hyperbolae in the second and fourth quadrants 
(assuming an energy above threshold), which is not illustrated here. 

At higher order these canonical coordinates are defined so that the Hamiltonian 
depends on the reaction coordinates through the reaction action I only. To estabhsh 
notation, we will write 

HiQ, P, q,p) = i7rc(/) + HUq,P, I) (2) 

where 

Hts{q,P, I) = Ho{q,p) + IH^{q,p) + PH^iq.p) + • • • 

The leading parts of Hrc{I) and i?o(?,p) contain the quadratic truncation of the 
Hamiltonian shown above. That is 

H,^(I) = -A7 + h.o.t. 

and 

d 

HQ{q,p) = + +h.o.t. 

i=l 

We refer to H^c{I) and Hts{q,p,I) respectively as the reaction coordinate part and the 
transition-state part of the Hamiltonian, whence the subscripts. 

We remark that for the following calculation to work, it is not necessary to put 
the transition-state part of the Hamiltonian in normal form. That is, the functions 
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Hm{.1iP) need not be written as functions of the transverse actions Jj = {qf+pf)/2. This 
allows greater latitude in the normal form construction and might in principle alleviate 
problems with small denominators. We also remark that since we will be looking later at 
complex solutions to the equations of motion, there is an implicit assumption throughout 
this paper that the Hamiltonian is analytic in the transition state region. 

2.2. The quantum normal form 

The transformation to the classical normal form described above is achieved by making 
a canonical change of coordinates so that the Hamiltonian takes the desired form. The 
corresponding procedure in quantum mechanics is to use a unitary change of basis to 
the same effect. While a canonical transformation does not uniquely define a unitary 
operator in the quantum formalism, it is a well-established feature of the quantum- 
classical correspondence that such a connection can be achieved within semiclassical 
approximation. The connection was outlined by Miller in using generating functions 
to write explicit approximations for corresponding unitary operators up to corrections of 
relative order 0{h). Explicit and concrete rules describing how unitary transformations 
may be constructed which achieve normal form to higher order in h have recently been 
published by Cargo et al in ^l] and used there to derive compact higher-order Bohr- 
Sommerfeld rules. In the context of critical transmission, which is the application of 
interest here, transformation to quantum normal forms have been exploited by Colin 
de Verdiere and Parisse in ^21 1^ provide connection formulas describing the 

behaviour of wavefunctions near hyperbolic fixed points, and these have been used in 
[TH] to calculate multidimensional quantisation rules which are valid near degenerate 
tori. 

An alternative approach is to work directly with the quantum problem in 
making the transformation to normal form rather than simply quantising the classical 
normal form as we do here (see [ini El CHI, for example). We do not adopt this 
viewpoint because we will in any case later need to employ semiclassically constructed 
nonperturbative unitary transformations to connect the normal form basis to the 
asymptotic basis used for the scattering matrix and there is no overall advantage in 
avoiding their use at this stage. 

In this work we are interested only in constructing the quantum normal form to 
leading order semiclassically — that is, neglecting corrections of relative order 0{h) 
in wavefunctions or terms of order 0{fi^) in classical symbols. This is achieved using 
the "preliminary transformation" in the language of ^I] and the resulting leading- 
order quantum normal form can be written straightforwardly as a direct copy of the 
classical normal form. A detailed discussion of this point would unnecessarily elongate 
the presentation here and will simply assert a direct equivalence (modulo higher order 
corrections) between quantum and classical Hamiltonians whenever necessary, referring 
the reader to jTH [T2] for a proper explanation. 
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We therefore start, in analogy with with a Hamiltonian of the form 

H = H,,{f)+H,,{q,pJ) (3) 

where / and Hts{q,p,I) respectively denote quantisations of the classical symbols / 
and Hts{q,p, i). There are ordering issues in this correspondence, of course, but for 
the purposes of making semiclassical approximation to leading order in h, it suffices to 
let / and Hts{q,p,i) be Weyl quantisations. The key feature here is that / and {q,p) 
act on different degrees of freedom and therefore commute. For concreteness, it may 
occasionally help to suppose that the transition-state part can be expanded in the form 

HUq,P, i) = Ho{q,p) + iH,{q,p) + PH^{q,p) + ■■■ (4) 

where we may in particular assume that 

[/,i/^(g,p)] = 0. 

The central result in this paper will be stated in an invariant way that does not refer 
explicitly to the normal form construction and the details of how this transformation is 
performed are not needed to use it. In addition, the essential idea of the calculation is 
understood simply on the basis of the normal form Hamiltonian itself. We will therefore 
simply quote the result and refer to Refs. jHH HH |T2I for detailed discussions of various 
approaches to making this transformation in practice. 

3. Scattering solutions of the quantum normal form 

Because higher-order terms in ^ couple the reaction degree of freedom to the transverse 
degrees of freedom, the normal form Hamiltonian is not separable in the simple-minded 
sense of the eigenvalue equation separating into a function of the reaction coordinate 
plus a function of the transition-state coordinates. Since the Hamiltonian depends on 
the reaction coordinate only through the reaction action /, however, it turns out that 
the eigenvalue equation nevertheless admits solutions which have a separable structure. 
We will use this property to reduce the transmission problem to one that is effectively 
one-dimensional and therefore solvable by standard techniques. Technical details of this 
reduction are given in the present section. In the next, it is shown how the results can 
be formulated in such a way that they no longer rely on an explicit consideration of the 
normal form. 

3.1. The reaction coordinate part 

The normal form construction provides us with a transformation to coordinates (Q, P) 
such that I takes the form 

i = \{QP + PQ). (5) 

which is the Weyl quantisation of QP. For interpretation of the results below in terms of 
conventional calculations, it may help to suppose that the reaction coordinate part H^-c 
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of the total Hamiltonian acts on functions of a coordinate x (with conjugate momentum 
Px) so that 

i = I{x,p,) (6) 

and that the problem in the ^-representation is close to a standard barrier-penetration 
problem. We could, for example, let {x,Px) coincide with the coordinates {qo,Po) defined 
in Section [2.11 as a rotation of {Q,P) in the reaction-coordinate phase plane. In that 
case the transformation from © to is achieved using a metaplectic rotation which 
rotates the phase plane clockwise through an angle 37r/4 (to give Figure EI). 

The operator I has continuous spectrum and in x-representation we write the 
(improper) eigensolutions in the form 

i^x{x) = X^j(x). 

Then 

Erc{I)=H,,{I = I) 

is the corresponding reaction-coordinate energy. These eigenf unctions are two-fold 
degenerate since we can send incoming waves from either the reactant or the product 
side of the barrier (see [Appendix A ). We will restrict our attention here to states which 



have an incoming component on the reactant side and outgoing components on both the 
reactant and product sides, but no incoming component on the product side, in which 
case there is a unique solution for each I. 

Either as an inverted parabolic barrier [H] or in the representation implied by (0) 
[T^ . the problem of finding eigensolutions of / can be solved exactly and solutions of 
the scattering problem written in closed form. Details are given in [Appendix A For 



present purposes it is sufficient to note that there is a simple relationship describing the 
relative fluxes in the incoming and outgoing channels. Let the scattering state ipii^) 
be normalised so that the incoming flux on the reactant side is normalised to unity (by 
construction, the incoming flux on the product side is zero). Then the outgoing fluxes 
on the reactant and product sides are, respectively. 

In the barrier-penetration picture, T(X) and -R(X) respectively represent probabilities 
of transmission and reflection. Note that by writing the transmission probability in the 
alternative form 

the unitarity condition 

R{X) + T(J) = 1 

becomes self-evident. 

We make the following observations concerning this result. 
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• There is a symmetry between -R(X) and T(T) on changing the sign of X. This is 
to be expected on the basis of a phase-space portrait in (Q, P) coordinates (see 
Figure Ej) in which changing the sign of X simply exchanges reactants for products 
in the outgoing channels, but is less obvious in a barrier-penetration picture. 

• The transmission and reflection coefficients do not change if we replace / by a 
Hamiltonian H{I) which is an arbitrary function of I. This will be obvious after 
generalised fluxes are deflned in the next section and a multidimensional version of 
this observation will be important in getting a simple formulation of the results in 
this paper. 

• The expressions in (0) give semiclassical approximations to transmission and 
reflection coefficients for a generic potential barrier and can be derived from 
the standard representation of the Schrodinger equation using the method of 
comparison equations f^. In the current calculation they are exact for any 
Hamiltonian which can be written as a function of the reaction action I alone. 
However, there is in general semiclassical error arising from the transformation to 
normal form in the first place, during which terms of 0{ff) arise in the Hamiltonian 
which are neglected in the current analysis. 

3.2. The transition state part 

In the transverse degrees of freedom corresponding to we suppose a discrete 

spectrum parametrised by X in the following way 

HUq,P,I)\Ml)) = E^,iX)\M^)). 

Here we suppose that a partial symbol Hts{q,p,X) is defined by replacing / by its 
eigenvalue X in Hts{q,p, /). If Hts{q,p, /) is given as a series of the form (jH) then we can 
write, concretely, 

Hts{q,p,X) = Ho{q,p) +XHi{q,p) +X^H2{q,p) + ■■■. 

In many interesting chemical applications, problems arise which have a Morse or Van 
der Waals type potential in the transverse degree of freedom for which the spectrum 
of Hts{q,p,X) is discrete at the bottom but becomes continuous above a threshold. 
We will confine ourselves, however, to energies at which asymptotically propagating 
scattering states correspond to reactant molecules in bound states, and for these 
cases the continuous spectrum (of Hts) does not participate. Rather than adjusting 
notation here to incorporate the continuous part of the spectrum we simply suppress it 
notationally and consider states labelled by the discrete index k only. 

Results like those in [8] can be obtained by ignoring the X-dependence of the 
states \ipk(X)) and approximating them by 1^9^(0)) (or equivalently keeping only the 
leading part HQ{q,p) in the expansion above), but here we want to investigate the effect 
of coupling between the reaction and transverse degrees of freedom seen in the full 
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Hamiltonian. For a fixed X we can assume that these states form an orthonormal set 
but note that we should assume in general that 

if X 7^ X'. This is the main point comphcating the following analysis and means that we 
should be wary of assuming "obvious" results when describing issues of normalisation. 

We will now describe how the discrete eigensolutions of the transverse problem 
combine with the scattering solutions found in the reaction-coordinate degree of freedom. 

3. 3. Eigenf unctions of the total Hamiltonian 

With the conventions described in Sections 13.11 and 13.21 

are eigenstates of the full Hamiltonian H in mixed position-bra-ket notation. These 
solutions satisfy 

H\^j,k)=Ek{I)\^i,k), 

where 

Note that the states |\l'j,fc) are separable in form even though the Hamiltonian itself 
is not strictly speaking separable, as discussed at the beginning of this section. This 
nonseparability manifests itself through the dependence of the transverse eigenstates 
|(y9fc(X)) on X and in particular through the nonorthogonality condition (|H)). 

Our aim is eventually to express results in such a way that explicit reference to the 
normal form transformation is unnecessary. For this purpose it is preferable to label 
states with the total energy instead of the reaction action, since the energy is defined 
independently of the representation used. For each value E of the total energy let, Tk{E) 
be defined implicitly solution of 

E = Ej,{l). 

Although we cannot write explicit expressions for Xk{E), we can suppose that in an 
energy range around threshold (where -Erc(X) = — AX+ 0(X^)), these functions are well 
defined and single-valued for each k. We then define scattering states labelled by the 
total energy E and the mode number k as follows. Let 

\^E,k) =i^i,iE)ix)\ME)), 

where for short we write 

Note that in view of (jHl) we have 

{ipkiE)\ME)) ^ 5kk' (9) 
since on changing k the action eigenvalue X = Ik{E) changes. 
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It is possible to normalise these states so that 

{<ilE',k'\'^E,k) = S{E-E')6kk' 

in the usual way, but this convention turns out not to be particularly useful for our 
purposes and we will not apply it. Instead we will normalise these states so that they 
have unit flux in the incoming reaction channel. A discussion of normalisation by flux 
will also be necessary to appreciate how arbitrary scattering states may be constructed 
from these separated solutions, so we will defer further consideration of such issues until 
a method of flux calculation has been outlined in the next section. 

4. Fluxes and sectional inner products 

Since the normal-form Hamiltonian is not of kinetic-plus-potential type, we cannot use 
the usual definition of current 
h 

J = (^*V^ - W^*) 

2im 

to calculate fluxes. There is a simple generalisation, however, which works for arbitrary 
Hamiltonians. 

Let B be a Hermitian operator which projects to one side of a section S which has 
codimension one in phase space. For example, if S is defined by fixing a configuration 
space coordinate then in position representation B can represent multiplication by the 
characteristic function of a region which has boundary S. More generally, B can be an 
operator for which a classical symbol such as the Weyl symbol rises from zero to unity 
in a classically small strip around S. Then the flux of a state |\E') across S is 

F={^//^^), (10) 

where we define a sectional overlap by 

{^//^^) = ^{m,Hm- (11) 

The notation here is adapted from JHI although a factor of ih has been introduced which 
will simplify matters later. Flux defined in this way is an integral part of transition- 
state calculations in the chemical literature (see [2II1 E] for example) and similar ideas 
are used in [TU It is easily verified that if B represents multiplication by the 

characteristic function of a region in configuration space then the flux reduces to the 
standard case of a surface integral of the current J over the boundary. Although a flux 
calculation requires only the diagonal case in (fTT|) . it is useful to allow the nondiagonal 
case in the definition of sectional overlap. We will find in particular that the space 
of scattering solutions with a given total energy E can be identified with a quantised 
surface of section in one of the channels and the sectional overlap then provides a natural 
inner product for the corresponding Hilbert space. 
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4.1. Fluxes for the separated scattering states 

We will now apply this generalised construction to compute fluxes in the normal form 
representation. Let the section E be chosen so that G) can be constructed in terms of 
the {Q, P) operators alone and commutes with q and p. In phase space this means 
that G) projects onto a region of phase space defined by a subset of the QP plane and 
independent of the {q,p) coordinates. Such a choice is natural if we view the complete 
system as a pinched waveguide (Figure in which a section obtained by fixing reaction 
coordinates is used to define fluxes in and out of the reactant and product channels. 




Figure 3. Sections wliich measure in and out fluxes in the product and reactant 
channels are illustrated schematically. Sections in full phase space are defined by 
fixing cither the coordinate Q or the coordinate P appropriately in each case. 

With the generalised definition of flux in (jlU|) we are not, however, confined to 
fluxes across surfaces in configuration space and are free to define sections in phase space 
which distinguish incoming from outgoing flux in the reactant and product channels. For 
example, referring to Figure a section defined by a vertical line in the right-half 
of the Q-P-plane measures incoming flux in the reactant channel, whereas a horizontal 
line S^"* in the upper half plane measures outgoing flux in the reactant channel and a 
horizontal line Ep"* in the lower half plane measures outgoing flux in the product channel. 
Fluxes across sections such as these corresponding to horizontal and vertical lines in the 
QP plane are especially easily computed in the QP representation. Details are given 
in [Appendix A| for the one-dimensional scattering solutions described in section 13.11 — 
here we simply note that relative fluxes are found to be of the form given in ((7j). Fluxes 
in the full system can be understood on the basis of these one-dimensional calculations 
in the following way. 

By an abuse of notation, let us denote a one- dimensional flux for one of these 
reaction-coordinate sections by 

{i^i/j:^Px) = -Ui^x\[QjMi). (12) 
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Although we use similar notation to (|l(Jj). it should be emphasised that this one- 
dimensional flux differs in having / and not a Hamiltonian in the commutator and 
is not therefore a straightforward physical flux. There is also a minus sign, which 
compensates for the fact that I{Q,P) generates a flow that is opposite in direction to 
the physical flow generated by the Hamiltonian — corresponding to E'f^{I) < in the 
discussion below. Although distinct from the physical flux, this quantity is prominent 
in the physical answer and the notation is useful for that reason. 

We will see that in the full space the separated scattering states \^i,k) have sectional 
overlaps of the form 

= -E'^{T){iPxU^i)hk'. (13) 
Alternatively, using states labelled by the total energy and normalised so that 

\i>E,k) = i^x,(E),.), (14) 

we have 

{^E,k' = {i^i, /M5kk'. (15) 

The main point here is that the multidimensional scattering states have fluxes, and 
therefore probabilities of reflection and transmission, which reduce exactly to the one- 
dimensional case and we can apply ((Zj) to scattering states of an arbitrary Hamiltonian 
in normal form. The transmission and reflection probabilities of the states \'^E,k) can 
therefore be written 

and we can now use this as a basis with which to treat reactivity of an arbitrary 
stationary state. 

We should emphasise that in view of the nonorthogonality of the transverse modes 
as expressed in (jSI) and Q, and in contrast to separable problems, these identities 
extending the one-dimensional results are not at all obvious. To prove them we consider 
separately the diagonal and nondiagonal cases, and treat them as follows. 

Derivation of (113^ in the diagonal case 

We treat the diagonal case k = k' first. Assume that can be expanded in the form 
(jH) and use 

[6, = [6, i]!""-' + i[e, + • • ■ + i""-^ [e, /] 

to deduce that 

{^T\[e,n\^i) = nii^-'{^i\[e,mx) 

and therefore that 

{^i,k\[^,Hts]\'fj,k) = J2mI"'-\ijj\[Qj]\ijj){M^)\Hr^{q,p)\Ml)) 
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= (^x|[e,/]|^x)^^, 

where in the last hne we have invoked the Feynman-Hellman theorem. On adding a 
similar calculation for H^c we get the claimed result. The important feature is that the 
factor -Efc(X) does not depend the choice of section across which to measure flux and 
relative fluxes reduce to the one- dimensional case. 



Derivation of (|13p in the nondiagonal case 

To treat the nondiagonal case we use the identity 

= + r + . . . (^^1 [6, /] 

■Tm 'Tim 

= (^xl[e,/]|^x) 

to deduce that 

-7"m TT -Tim TT 



X-T 



On summing over m and doing a similar calculation for the reaction-coordinate part of 
the Hamiltonian we find that 

In terms of the energy-labelled states \^E,k) = \^ik{E),k)i we have 

{^E>,k> //T^E,k) = ~ x^{E) -Xk'{E') ^'^^^'^^'^ /s'^TkiE)){fk'{E')\ifk{E)). 

These sectional overlaps vanish as E' approaches E unless k = k', in which case we 
recover the diagonal result. This completes the derivation of 

4-2. Fluxes for general scattering states and reduced Hilbert space 

Let denote the subspace of scattering states of a fixed total energy E which are 
incoming on the reactant side and have no incoming flux on the product side. In the 
applications we consider the scattering problem will asymptote to a waveguide-type 
problem in which there are a finite number M of propagating channels (corresponding 
to the energetically accessible bound states of the reacting molecules) and the space 
will therefore be finite-dimensional. We assert that the space is in fact essentially 
a quantisation of a classical surface of section (in the incoming reactant channel) in the 
sense of Bogomolny j2S] , with the inner product being given by sectional overlaps of the 
form in (jllj) . 
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Let us write a general asymptotically-propagating scattering state in the form 

M I T \ 

i$)=i:«^7===- (17) 

k=l \J mJ T,^rAl)x^) 

Then the total incoming flux for this state is 

M 



k=l 



while in view of ()15|) and ()16p the outgoing flux in the reactant channel is 

M 



M 

^ i 

k=l 

and the outgoing flux in the product channel is 

M 

^ i 



27rXfc(£;)/;il 

fc=l 

Furthermore, between any two such scattering states the sectional overlap 

M 
k=l 

provides a natural inner product for TH}^. 

We can regard Tt^ abstractly as a space spanned by an orthonormal basis {\k)}'^i^^ 
whose elements 

\^E,k) 



\k) 



are in one-to-one correspondence with the scattering states \^ E.k)-, normalised to have 
unit incoming flux. A general element 

M 

\v) = ^ak\k) 

k=l 

of this reduced space can be extended to a scattering state of the form given in (fT7|) for 
which the total energy is fixed but which is not necessarily an eigenstate of the transverse 
Hamiltonian Hts- The inner product between two reduced states can be defined through 
sectional overlaps 

of the corresponding extended states. 

The benefit of this abstraction is that we can compute outgoing fluxes in the product 
channel using matrix elements 

Fr = {miEM (18) 
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of a reaction operator 

^(^) - E 1 = E 1 +e-^.x.,.,/. l^->(^i (19) 

fe=l fc=l 

which is diagonal in this basis. By writing the outgoing flux as a matrix element of an 
operator defined on IH}^ we have in large part achieved the goal of this paper, which 
is to generalise the construction in so that there is no longer a restriction to states 
supported in a classically small region of phase space. In fact, the only restriction on 
incoming states here is that they should be supported in the region of phase space 
where the normal form in (0) provides an accurate description of dynamics. Although 
undoubtedly a restricted subset of phase space, this domain has classical dimensions. 

The current version is tied to the normal form, however and in order for this operator 
to be of any practical use, we really need a way of constructing it which does not call 
on us explicitly to construct the normal form or the basis vectors \k). As a start in this 
direction, note that if we denote 

f{E) = ^e-2-^^(^)/^|A;)(A;|, (20) 

k 

then TZ{E) can be written in the form promised in the introduction. This is precisely 
the form given in 0, where T^E) was a tunnelling operator defmed as a quantised surface 
of section map in the neighbourhood of a complex periodic orbit. We will show in the 
Sectional that the same interpretation can be imposed on 'T{E) in the present case with 
the difference that, unlike in jH], restriction to a neighbourhood of the periodic orbit 
small enough for linearised dynamics to be used is no longer necessary. 



4-.3. The microcanonical cumulative reaction probability 
The trace 

M 

N{E)=TrnE) = j: (21) 
k=i 

of the reaction operator is the so-called microcanonical cumulative reaction probability. 
Results for N{E) equivalent to those that would be obtained by using linearised 
dynamics in 'R-{E) in the manner of 8j were obtained by Miller in ||24j. Semiclassical 
approximations for N{E) have also been given in [201122] which include nonlinear effects 
in the transition state degrees of freedom by using expansions in the transition-state 
quantum numbers (a related treatment of tunnelling using normal-form coordinates 
has been given in [2n])- A thermalised version has been given in |2Z] and see |23 fo'^ 
a semiquantum calculation. A discussion emphasising the fluctuations that occur in 
N[E) as the summands in (f?T|) switch on with increasing E can be found in l3Uj. 
We also remark that a discussion of N{E) has recently been given in JB] which uses the 
same language of normal forms that we use here. 

The major benefit of the current work is that, once a prescription has been given in 
the next section for computing 'R-{E) without recourse to the normal form, we will have 
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an explicit prescription for distributing the total reaction probability in N{E) among 
incoming states or, equivalently, using the Wigner-Weyl formalism, for understanding 
how the reaction probability is distributed in phase space. In principle, this calculation 
can be implemented simply by computing complex trajectories near the complex periodic 
orbit used in [Hj and does not require a particular deconstruction of the Hamiltonian 
such as provided by a normal form (although normal forms do help in inverting the 
operator 1 + T{E) in ((H) as we will discuss). 

5. Getting away from the normal form 

Equation ((T)) promises the ability to treat reaction problems without having to deal 
explicitly with the normal form and separated stationary states. In this section we show 
how this can be done by interpreting 'T{E) as a tunnelling operator in the sense of jHllSl 
which can be understood independently of the normal form construction. 

5.1. The complex return map in normal form coordinates 

We first describe how a complex return map is expressed in terms of normal form 
coordinates. Let (/, 6) be action angle variables in the QP plane such that in the 
positive quadrant 

P = Vle-r 

We have / = QP, consistent with previous notation. Although the system is hyperbolic 
and not normally associated with periodic motion, there is in fact an imaginary period, 
expressed by the identities 

Q{e + 2'Ki) = Q{e) 
p{e + 2'Ki) = p{e). 

In other words, the Hamiltonian flow generated by / in the complexified QP plane has 
period 2Tri. We will now show that an extension of this periodic flow to the full system 
can be used to define a complex Poincare return map and that this map has a fixed 
point corresponding to a complex periodic orbit. 

Let us restrict initial conditions to a surface of section defined by fixing Q as 
in Figure El along with the total energy E, and let the condition 

E = H{q,p,I) (22) 

implicitly define the function h{q,p,E) on by 

I = hiq,p,E)-eiE). 

Here we regard {q,p) as canonical coordinates for S"^ and e{E) is defined so that the 
minimum of h{q,p, E) on S"^ is zero, as described more explicitly below. Differentiating 
(|22|) with respect to the transverse coordinates q and p while keeping E fixed gives 
f)M 

= -T^Vh + VH, 
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which in turn gives 

/ dg \ 
- "^^^ - " dp ' 

where Xh denotes the projection of the full flow vector Xfj defined by H onto the 
transverse degrees of freedom {q,p) and we have used 6 = dH/dl. The flow defined by 
Xh can therefore be regarded as a restriction to the {q,p) degrees of freedom of the full 
flow, reparametrised so that time t is replaced with the angle variable 6. 

Letting 6 evolve from to the flnal value 27ri, trajectories are described in full phase 
space which start on S'^ and return to it — recall that the [Q, P) coordinates which are 
used to deflne S'g are periodic under this evolution. Integrating the flow vector —Xh 
for a time 27ri then generates a complex symplectic map 

-r yin v-iin 
• ~^ ^R 

which we can denote by 

= exp[-27TiXh] (23) 

in Lie-algebraic notation. This is precisely the classical map used to construct the 
tunnelling operator in |H]. 

Before describing explicitly how the quantisation works, it is helpful to see how 
the complex periodic orbit which provides a flxed point of JF arises in normal form 
coordinates. By construction, the quadratic part of the Hamiltonian H{q,p, I) is elliptic 
in the transverse degrees of freedom. As a result, the sectional Hamiltonian h{q,p,E) 
has a minimum {qf.{I) , Pe{I)) , for sufficiently small / at least, for which 

Vh{q,{I),Pe{I),E) = 

and which coincides with the origin of S'^ in the threshold case / = 0. We define e{E) 
above so that h{qe,pe,E) = and a Taylor expansion of h{q,p,E) about {qe{I) , Pe{I)) 
begins with quadratic terms. This minimum is an equilibrium of the fiow defined by Xh 
on S"^ and is therefore a fixed point of J-". In full phase space, the trajectory starting with 
coordinates {qe{I),Pe{I)) on S'^ evolves so that the coordinates / and {q,p) are fixed 
and defines a periodic orbit as 9 evolves from to 2Tii. The time period corresponding 
to this evolution is —ir{E) where 

T = = -T^^ TTT^ = 2'ne'iE). 

9 dH{q,,p,J)/dI ^ ' 

and its action is an imaginary number S{E) = iKQ{E) where 

Ko = - ^ Id9 = 2nl = -2ne{E). 

In jH] a linearisation of dynamics about this complex periodic orbit was used to 
approximate 'R-{E) which here corresponds to truncating a Taylor series of h{q,p,E) 
about {qe{I) , Pe{I)) at quadratic order. The essential conclusion of this paper is that a 
complete description of the reaction operator can be achieved simply by replacing this 
truncation with the full sectional Hamiltonian h{q,p,E). 
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5.2. The tunnelling operator 

The tunnelling operator is defined to be a quantisation of the classical map J-' jHl IT9] . 
Using we can write concretely, 

where h — e is the restriction of the operator / to the quantum analogue TiJ^ of SJ^ 
defined in the Section EJ In that case 

Ch-e)\k)=Ikmk) 

and 

and the identification in the previous section of T{E) in (pUj) as a tunnelling operator 
is confirmed. 

It should be emphasised that while the sectional Hamiltonian h was used in 
making this identification, it is not necessary to construct it explicitly, or even to 
refer to it, in order to construct the map and to approximate its quantum analogue 
'f[E) semiclassically. The map JF can be constructed simply by integrating orbits as 
described in references jSl E] ■ From the dynamical characteristics of these orbits. Van 
Vleck-type approximations for 'T{E) can be written as described in [22], for example. 
Alternatively, the Weyl symbol of T{E) can be obtained as described in jHIj. Coherent 
state representations are also possible ([22123 and references therein). In all of these 
approaches the djTiamical information needed is naturally provided as a result of the 
orbit computation and h{q,p) is not needed explicitly. 

We note finally that while it would be quite easy to write a version of p9|) that 
describes the full scattering matrix in the normal form representation, so that we could 
determine the distribution of product states for each reactant state, it is less obvious 
how the normal form result could be interpreted in a basis-independent way in that 
case. Any such reformulation would have to take into account the fact that orbits 
contributing to the scattering matrix itself ^0] depend singularly on initial conditions 
near the reacting boundary. The simplicity of the normal form representation suggests, 
however that such a uniformisation might be feasible, although we do not pursue it here. 

5.3. Asymptotic basis for the reaction operator 

Although the reaction operator is diagonal with respect to the stationary states \'^E,k) 
computed in terms of the normal form, it will not in general be diagonal in the basis of 
asymptotically decoupled stationary states that is used to write the scattering matrix. 
In fact, since the normal form will in general only provide an accurate description of 
dynamics in a neighbourhood of the transition state, we need an independent method 
to describe how these states can be extended to the asymptotic regions of the reactant 
channel and the operator 'R-{E) written in the standard asymptotic basis. 
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To achieve this we note simply that once the transmission problem has been solved 
locally in the transition state region, the solution can extended to the asymptotic region 
by applying quantised surface-of-section maps, such as described in j^, for example. 
Once outside the transition state region, these mappings can be constructed on the basis 
of primitive semiclassical approximations and amount in the classical picture simply to 
conjugating the complex Poincare map J-' with standard real ones. As when transforming 
to the quantum normal form in the first place, the end result of this process is easily 
stated and the details omitted in the interests of brevity since they follow discussions 
elsewhere 

To be more specific about this conjugation, let the coordinates {x,y,px,Py) be as 
described in the introduction and let the Hamiltonian decouple asymptotically in the 
reactant channel according to 

H{x,y,Px,Py) ~ H^y^p{y,py) + Htr{x,px). (24) 

Here Hasyrap{y,Py) describes the internal vibrational motion of the reacting molecules 
and Ht^{x,Px) is a kinetic energy term for the relative motion of centres of mass. In the 
simplest atom-diatom collinear case we might have 



2M 

where M is the atom-diatom relative mass and 

H^jmp{y,Py) = ^ + ^ABiy) 

where m is the diatomic relative mass and Vab{(i) the diatomic interaction potential. 
The separated asymptotic solutions are denoted 

\^n,E) = XE,n{^)\i'n) 

where XE,n{x) is a plane- wave eigenfunction of Htr with energy E — En and normalised 
to have unit incoming flux while \ipn) is an eigenstate of internal dynamics with energy 
En- There is a phase convention implicit in writing the scattering matrix, which in the 
present case amounts to specifying the phase of XE,n{x). In the standard plane- wave 
case 

XE,n{^) = 

\/Vn 

where f„ = hkn ensures unit incoming flux, we can understand the phase of Xs,n(a^) as 
being fixed at x = 0. 

The phase convention in the asymptotic regime is then obtained by extending the 
complete scattering state \^n,E) from a section S^j defined by a; = to a section Xl^-^™^ 
defined by a large and negative value of x. This extension is affected in the semiclassical 
scheme by quantising a surface of section mapping 

constructed using the uncoupled dynamics of the Hamiltonian in ()24|) — since the plane 
wave is travelling to the right, this will be achieved in negative time. Scattering of the 
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f lm(t) 



Figure 4. The path in the complex time plane which achieves the conjugation H25II 
is illustrated schematically. The dashed lines indicate evolution under decoupled 
dynamics. For the complex periodic orbit defining a fixed point of the contour 
segments corresponding to JFqj -^asymp and their inverses are parallel to the real axis, 
while the segment for J-'ts is parallel to the imaginary axis. The time evolution of each 
segment develops both real and imaginary parts as the initial condition moves away 
from this fixed point. 



resulting incoming state then proceeds by first mapping from the section j]^^™^ to a 
section near the transition state using the fully coupled dynamics, 

.77 . yasymp yts 

-'"asymp • ~^ ^R- 

If S^^ is in the domain of the normal form then its transmission probability is understood 
using a complex return map 

y-its y-its 

•r ts • — ^ ^R 

of the kind described in Section I5.ll Finally, we complete the transformation to the 
asymptotic basis by mapping back to the asymptotic section j]'^^^^^ using the inverse 
•^asymp -^asymp and then mapping to S^j using the inverse J-'q^- The end result is a 
map 

T = Tq ^asymp-^ts^asymp-^O (25) 

which is conjugate to jF^g but which is adapted to the phase convention used for 
the scattering matrix. A quantisation of this map gives a tunnelling operator TiE) 
appropriate to the asymptotic basis defined by the states | 

Note that by conjugating J^ts by J^asymp, we can extend the return map far beyond 
the domain where the normal form applies. The outer conjugation by JFq is not strictly 
necessary to calculate reaction probabilities of incoming states | since the transverse 
parts I'l/'n) are eigenfunctions of the quantisation of JFq and the resulting eigenphases 
cancel in the transformed version 
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of (jl8|) . However phases are important for cross terms if we want to treat general 
scattering states 

M 
n=l 

and they are also important for representations of 'R-{E) in phase space. We also remark 
that J^^isymp will not in general commute with Tts and the reaction operator will therefore 
not be diagonal in the asymptotic basis defined by the states On noting that 

the conjugation in (^3)) amounts simply to a change of representation so that S^j can 
be identified with the section E'g, however, we see that the tunnelling operator written 
in this basis is not fundamentally different from the one described in Section 13721 

Finally, we note that the conjugation in is particular to a waveguide problem 
in which the plane wave part is written in terms of a coordinate x with phases fixed at 
X = 0. If different conventions are used for the asymptotic coordinates or for the part 
XE,n{,x) of the scattering state, then the part JFg of the conjugation must be redefined 
accordingly. 

5.4- Limitations of the derivation 

No approximations are made in getting to the reaction-operator form of the outgoing 
flux in (|T8|l and ()19|1 once the quantum normal form of the Hamiltonian in is written 
down. The sources of error are in transforming to the quantum normal form in the first 
place and in transforming to the asymptotic scattering basis after the reaction operator 
has been found in the normal form representation ()18|1 and (fTT^ . We now comment on 
some of the issues affecting this approach. 

The classical normal form is a formal series which describes the dynamics locally 
in a neighbourhood of the transition state region. It is clear that writing a quantum 
version of it as we do here is a formal step which will ultimately require a more careful 
justification. Issues of convergence become especially important if we pursue the limit 
^ in its literal sense and we have not addressed such questions here. From a purely 
practical point of view, however, if one's aim is to achieve an approximation that works 
well for a small but fixed value of then it suffices to describe the dynamics to a 
corresponding level of accuracy in phase space and the normal form is certainly capable 
of that in the sorts of parameter regimes that arise in chemical applications jSl]. 

Results of a numerical investigation are outlined in [33] which indicate that the 
expression in works very well when T{E) is computed directly from the complex 
orbits of the Poincare return map JF. We also note that since ((T)) can be interpreted 
theoretically without reference to the normal form, it seems natural to expect that it 
applies independently of the normal form itself. We therefore conjecture that, despite 
the limitations of the derivation presented here. Equation is in fact "classically 
exact" in the sense that no errors arise from classical dynamics side of the calculation 
once T^E) is interpreted as the quantisation of and the only approximation is the 
usual semiclassical one which vanishes as /i — > 0. 
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It is important to add the qualification, however, that even if the conjecture is 
correct, there are good reasons to expect it to apply only locally, at least in the 
simple form described in sections 15.11 and 15.21 The return map describes a unique 
image for initial conditions in a neighbourhood of the complex periodic orbit and for 
energies sufficiently close to threshold. In practical terms, this means there is a unique 
complex solution satisfying the boundary conditions required of orbits by semiclassical 
approximation of in the usual representations [SI Sufficiently far away, 

however, bifurcations are likely to occur where this structure breaks down and these 
are not described by the current formulation. In fact, it has been shown in (and 
see jSni inn IHHl inn] for related work) that complex orbits contributing to the scattering 
matrix can be chaotic and are subject to intricate pruning by the Stokes' phenomenon, 
while numerical evidence jH^ suggests that the same is true of the orbits contributing 
to TiyE) sufficiently far from the centre of the reacting region. At an even more basic 
level, the NHIM itself may undergo bifurcation once the energy rises far enough above 
threshold and in this case the whole bottleneck picture at the basis of our calculation 
is no longer correct. Global recrossing may occur and resonances arise in the quantum 
mechanics which are not described by the simple picture of transmission probability we 
have here. It should be stressed, however, that whatever the limits are on the domain 
where contributing dynamics are simple, they are independent of h and therefore have 
classical scales. 

Finally, we remark that even though the tunnelling operator T^E) can be routinely 
approximated using semiclassical approximations, the inverse of 1 + T(£') that occurs in 
is more problematic. Closed form analytic approximations are possible if we know 
the sectional Hamiltonian h{q,p) (see jH for the harmonic case) but more work is needed 
to provide an approximation that works directly in terms of the map JF. 

6. Conclusion 

We have characterised the semiclassical transmission of waves across a phase-space 
bottleneck using a reaction operator constructed from a complex Poincare mapping. 
In contrast to previous work [H], this construction is not restricted to energies and parts 
of phase space in a classically small neighbourhood of the transition state at threshold. 

A phase-space representation of the reaction operator will be largely supported in 
the classically reacting subset of phase space, but will also incorporate tunnelling and 
other quantum effects at the boundary of this region, where trajectories approach the 
NHIM along its stable manifold. The only dynamical information needed to apply the 
approximation described in this paper is contained in the complex Poincare mapping 
and explicit consideration of normal forms, or other special assumptions regarding 
the dynamics such separability or adiabatic approximation, are unnecessary. We also 
note that the particular trajectories used to define this map are well behaved at the 
reacting boundary. Therefore, despite the fact that the fate of trajectories changes 
discontinuously across the reacting boundary, the result here uniformly describes the 
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transition from classically allowed transmission inside the reacting region to reaction 
entirely by tunnelling outside it. 

We conclude by noting that in its current form the result here is restricted to 
collinear problems. In order for the approach to be used in completely realistic models, 
the marginally stable degrees of freedom associated with rotational symmetry will need 
to be incorporated. This aspect needs further investigation. 
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Appendix A. Scattering for the one-dimensional normal form 

The one- dimensional problem is especially easily solved in a representation in which the 
reaction action operator takes the form © and P acts on functions of Q according to 

In this case the eigenvalue equation lip = Tip is a first order differential equation 

and this can be solved in elementary terms without the complication of parabolic 
cylinder functions that arise in the conventional representation of an inverted oscillator 
jH]. This approach has in particular been exploited in [T2] to treat scattering in one- 
dimensional networks of tori and we refer to that publication for more detail of the 
following calculation. It is useful, however, to reiterate some of the main points here 
and to emphasise flux calculations, which form a basis for the discussion in the main 
text. 

The eigenvalue equation ()A.1|) leads to doubly degenerate eigenf unctions ipjiQ) of 
the forms 

and 

^iiQ) = ^U-Q) 

respectively. We will concentrate on the solution ip^lQ), which as we will now show 
represents bombardment of the equilibrium from the reactant side and is therefore the 
solution singled out in section ITTl 

The incoming flux is normalised as follows. Let the projection operator have the 
Q-r epresent at ion 

QHQ) = <diQo - Q)HQ) 
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so that it measures flux from right to left in the QP-plane across a section defined 
by Q = Qo- Then 

^[0, /] = ^ P] + [e, P]q) = -Q5(Q - go) 

and sectional overlaps are of the form 

ii^/^J) = -^(^|[e,/]|V) = Qol^(Qo)r. 

The case Qq > corresponds to incoming fiux on the reactant side. We then denote 
= E^^ and get 

n 

A section with Qq < gives a flux in the incoming product channel and this vanishes 
for the state ipj{Q), consistent with our interpretation of it as an incoming state in the 
reactant channel. 

Outgoing fluxes are naturally measured in momentum representation 
using projections of the form 
In this representation we have 

^[0, /] = ^ ([0, Q]P + P[^, Q]) = -PS(P - Po) 

and therefore 

To complete the calculation we therefore need to evaluate 

1 / I PI \ poo 



V27rh 



_g-io-7r/4+(T7rJ/(2?i)p / _ _|_ 11 j \p\-l/2-iI/h 



2^ V2 h 

where a is the sign of P. Note that <Pj{P) ~ const x |p|-i/2-j2:/'i j^^g ^j^g same 
dependence on its argument as found in the Q-representation, which is to be expected 
since there is a symmetry between Q and P in J. 

In either case, for flux calculations it suffices to know that 

awJ/h 1 

I + / D\ 1 2 

\'^l[^>\ = g7rX/R_^g-7rX//i, |p|' 

where we have used 

2 

TT 



cosh ttT/ h 



Semiclassical transmission across transition states 



26 



The J-flux across a section S"*^* defined hj P = Pq is therefore 

where here a is the sign of Pq. This is positive in the reactants-out channel (P > 0) 
and negative in the products-out channel [P < 0) which, when we remember that is 
defined so that upward fluxes are positive, is consistent with Figure |21 We have therefore 
confirmed Equation ((7|). 
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